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Abstract 

We consider the "second-class current" decay t~ — ► -K~nv T from several points of view. We first 
focus on the decay rate as expected within standard weak interaction and QCD due to isospin 
violation. The decay contributions divide into P- and 5-wave parts. The former can be reliably 
estimated using the pnir coupling inferred from the rates and Dalitz-plot distributions of n — * 3ir 
decays. The somewhat larger 5-wave part, which was previously computed using chiral pertur- 
bation theory, is estimated from a simple qq model. Both estimates of the 5-wave part depend 
on whether the ao(980) scalar particle is a qq or some other (4-quark) state. Finally, we discuss 
genuinely new, non-F — A scalar weak interactions. The t~ — ► ir~rju T decay provides information 
on this question, which nicely complements that from precision [3 decay experiments. In summary, 
we discuss the possible implications of putative values of the branching fraction B(t~ — > ir^nvr). 
In the case of larger values, in particular of the S'-wave part, not only will detection of the decay 
be more likely and more reliable, its implications will be more far-reaching and interesting. 



I. INTRODUCTION 



The weak decay t~ — > n~rii> r , an example of "second class current" decays introduced 
by Weinberg [1], may soon be observed or tightly bounded by the B factories. This isospin- 
and G parity-violating decay is suppressed by the small value of (m d — m u ) / Aq C d or a £M . 
Various estimates [2] using chiral perturbation theory or other methods have predicted this 
decay's branching fraction to be 

B = B(r~ -> n'TjUr) = (1.3 ± 0.2) x 10" 5 , (1) 

far below the present CLEO upper bound of 1.4 x 10~ 4 [3]. In view of the possibility of new 
measurements, we point out interesting consequences of various B values. 

The plan of this paper is as follows: In Sec. II we present the kinematics and some other 
general aspects of the r~ — > ix~r\v T . The contribution of the vector (L — 1) n^r] final state to 
B is discussed in Sec. Ill, assuming that the L — 1 and 1=1, n'r] final state is dominated by 
the p~ meson. Sec. IV addresses the contribution of the J p = + n~r] state to B. The analog 
of the p~ here is the 1=1, % (980) state, whose coupling to the vector current relates to a 
longstanding question on whether the a (980) is a ud state or a udss / K K-threshold state. 
The ud assumption was implicitly made in the chiral-Langrangian calculations predicting 
Eq. (1), where (980) dominance was used to analytically continue the calculation of low- 
energy decays to the r~ decay of interest. We briefly discuss another naive quark- model- 
based estimate. Sec. V addresses the possible relation between B and precise measurements 
of /?-decay spectra from trapped radioactive ions. Such measurements can be used to search 
for scalar interactions, in addition to the standard electroweak (V — A) ■ (V — A) interaction. 
In the concluding Sec. VI we present putative B values and/or bounds on B with implications 
for the discussions in the former sections. 

II. KINEMATICS OF THE r~ -*■ vr"^ DECAY 

Only the vector weak current Vjyx) = u(x)'j fM d(x) contributes to the hadronic part 
(0| J^\rin~^ = of the current-current interaction, since the 1 + and 0~ parts of the axial 
current cannot create natural-parity states of two pseudoscalars. The matrix element 
can be decomposed into a J p = + part and a 1~ part in the rest frame of the r]7r~ system 
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as follows: 

(0\V,\n V ) = f 1 (s)q IM + f (s)Q„ (2) 
where fi is the coefficient of the state with angular-momentum L, 

Q = q n + q v , 
q = a(s)q w - q v , 

S EE Q\ (3) 

q x is the four-momentum of particle x, and 

a s = " r: 4 

is chosen so that Q • q — 0. In the rest frame of the fpx - system, g is a space-like vector: 



q = (0, \q\ cos#, |g| sin#, 0), (5) 

where is the angle in this frame between q and the recoiling neutrino momentum. The 
L = and L — 1 amplitudes interfere in the angular dependence dT/d(cos(9), but not in 
the total decay rate obtained by integrating over d(cos(9), namely, 

^ = K 1 \f 1 (s)\ 2 + K \f (s)\ 2 , (6) 

with the K L being kinematic factors. Thus, either the S- or P-wave contribution yields a 
lower bound on the total rate. We proceed with an estimation of the magnitudes of these 
contributions. 



III. ESTIMATING THE L = 1 CONTRIBUTION OF THE tt"^ STATE 

The decay r~ — > 7r~7r°z/ T comprises 25.5% of all r~ decays, and is completely dominated 
by p~ exchange. Similarly, our estimate of the L — 1 contribution to the decay t~ — > ix~r\v T 
assumes p _ dominance, taking place via r~ — > p~v T followed by p~ — > rpT~. We thus expect 
the L = 1 component of £> to be 

\QpiTK J \Pp— >7T7T 

where ^ p7r and g p7r7r are the p — > ?pT and p — > 7T7t coupling constants, respectively, and the 
cubed ratio between the daughter momenta in the two decays is (pp^^/pp^^) 3 = 0.07. 
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B(r ^p i/ T ), 



(7) 



Since the decay p~ — > ^7r~ has not been observed, we obtain the coupling constant g vpn 
from the Dalitz-plot distribution of the decay rj — > 7r + 7r~7r° and the branching fraction 
B(rj — > 7T 7r 7r°). The three-pion Dalitz plot is customarily described with the variables 

y S |r„-i, (8) 

where T c is the kinetic energy of the pion with charge c, and 

Q = — 2m % + — m n o ps — 3m,,-. (9) 

Henceforth, we ignore the difference between the charged and neutral pion masses. The 
matrix element for r] — > 7r + 7r~7r° is taken to be the sum of a scalar and a vector exchange 
contribution, the latter dominated by the p(770): 

M+-o = M S + M P+ +M P -. (10) 

A p° contribution is forbidden due to charge conjugation conservation. Properly accounting 
for the number of diagrams and identical particles, the i] — > 7r 7r°7r° matrix element is 

3 

M 00 o = ~^ Ms - ( n ) 
The branching fraction of this decay gives the absolute value of the scalar matrix element, 

\M S \ 2 = 8(2n) 3 m v T v B( V - nW)^-^ = 0.065, (12) 

where we used the measured values of the rj mass, width, and 7r°7r°7r branching fraction [5] , 
the phase-space differential is dE 1 dE 2 = (Q 2 /6V3)dX dY, and Si = 2.75 is the area of 
the Dalitz plot. The scalar particle exchanged is assumed to be very broad, so that the 
distribution of events over the relatively small Dalitz plot is essentially uniform. 
We take the vector matrix element to be 

(P, + P ± ) • (P T - P ) 



A'lpT 9r)pn9pirTr 



(P T + P ) 2 -ml-iT 



2m v (E -E T ) 
9rm9p ™2m v E ± + M 2 -lmX 1 J 



where E + , E_, and E are the rj- rest-frame energies of the 7r + , tt , and n°, respectively, and 

M o = m l ~ \ m l ~ m l + lT P m P- ( 14 ) 



Replacing the energies with the Dalitz-plot quantities of Eqs. (8) and (9), the sum of the p + 
and p~ contributions is 



rY - \r 2 {Y 2 + X 2 ) 
M P -+M P+ = -2^^^ x2 



where 



, (- 

3 1 3 V3 

2 3 

rY + — (Y 2 - X 2 ) + — (X 2 Y - Y 3 ) 
o 9 



r = 



M 2 



0.14 + 0.03i. 



(15) 



(16) 



and the last line of Eq. (15) is obtained from a Taylor expansion to order r 3 . 

Squaring the sum of the scalar and vector terms, again keeping terms to order r 3 , we 
obtain 



where 



\M+- \' 
\M S \ 2 



1 + aY + j3Y 2 + 7 X 2 + 5Y 3 - 5YX 2 , 



(17) 



a = -4:g vp7T g pn7T ?R:{M* s r} 



12 ' 
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-^gr,p7v9pw7r^ {M* s r 2 } + A(g vp7r g p7Tn ) 2 \r\ 

^9r, P n9pnn^ {M* S T 2 } JJJ^: 



\M S r 



Q9vpir9fm*$l {M* s r 3 } + -(g^g^f'R (r(r 2 )*} 



\M S \ 2 ' 



(18) 



The product of coupling constants g^g^ is obtained by comparing the coefficients of 
Eq. (17) with the Dalitz-plot distribution of the decay rj — > tt + 7t~7t°. A high-statistics study 
of this distribution has been recently performed by the KLOE collaboration [4] , yielding the 
parameterization 



\M+^ \ 2 oc 1 - 1.09Y + 0.124F 2 + 0.057X 2 + 0.14F 3 . 



(19) 



We ignore the measured coefficient errors, as they are much smaller than the theoretical 
errors associated with our model. From the coefficient of the Y term in Eq. (19) and the 
first of Eqs. (18), one obtains the product of coupling constants 

1.09 M s 



9r) P Tr9pmr 



4 3?{r} 



= 0.51, 



(20) 
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where M.s was taken to be real. The accuracy of the model may be judged from the values 
it obtains for the other coefficients: 



|M + _ | 2 oc 1 - 1.09Y + 0.27F 2 + 0.05X 2 + 0.03F 3 - 0.03FX 2 . (21) 

Allowing M.s to have a complex phase does not improve the agreement between Eqs. (19) 
and (21) significantly. A related cross-check is provided by the ratio of branching fractions 
B(r] -> n + n-n°)/B(r] -> 7r 7r 7r ) = 0.70. The value predicted by Eqs. (17) and (12) is 0.71 
when using the experimental coefficients of Eq. (19), and 0.76 using those of Eq. (21). 
Taking the matrix element for the decay p — > nn to be 

M p = gpmr e®(P + -P-y, (22) 
the coupling constant g p7T7T is determined to be 



9pn, = J 3 P P = 6.0. (23) 

y Pp— >7T7r 

Eqs. (23) and (20) then give 

g Vfm » 0.085. (24) 

A similar calculation by Ametller and Bramon [6] yielded the ratio g VP n/9pmr — 0.011±0.002, 
consistent with our results. 

From Eqs. (7), (23), and (24) we calculate the L — 1 component of the t~ — > ix~r\v T 
branching fraction, 

B L=1 w 3.6 x 10" 6 . (25) 

We also obtain 

B(p - rjn) = 9 l py ^ » 1.4 x 10" 5 , (26) 
67rm^r p 

far below the current experimental limit of 6 x 10~ 3 [5]. 



IV. THE L = CONTRIBUTION 



The Contribution of the (L = 0) n^r] state to B is not as readily accessible to a phe- 
nomenological estimate as that of the L — 1 state. The observed p~ dominance in the 
7r~7T° final state of the r~ decay is expected, since the p~ has the quantum numbers of the 
hadronic vector current u^ p d. It is therefore natural to assume that it also dominates the 
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(L — 1) n^r] final state, although this decay is suppressed by isospin violation. This is not 
so for the superficially analog case of Oq (980) and the scalar contribution to £>. In Ref. [2], 
the Oq (980) dominance of the (L = 0) n^i] channel in weak decays was used to extrapolate 
the low-energy amplitude for rj — > 7r~e + z/ e (computed via chiral perturbation theory) to the 
decay r~ — > ~K~r\v T and obtain the estimate of Eq. (1). The resulting scalar contribution to 
B is then ~ 3 times larger than the vector contribution. This extrapolation is questionable 
not only because of the large change in Q 2 from ~ 0.15 GeV 2 to ~ 1 GeV 2 . The key point 
is that a^(980) (just like its 7 = counterpart / (980)) may well be a four-quark udss 
state, a view suggested early on [7] and adopted recently by the Particle Data Group [5]. In 
this case, the ao(980) coupling to the ud scalar current is "Zweig-Rule" suppressed, and the 
four-quark state will not dominate the decay in question. 

Several considerations suggest that the ao(980) and /o(980) states have significant four- 
quark contributions: 

1. The widths r(/o(980) — > titi) ~ T(ao(980) — > ni]) ~ 50 MeV are anomalously small 
for an S-wave qq state. Since the lighter, 770-MeV p has a P-wave decay width of 
150 MeV, the a (980) / (980) and widths should have been vastly larger. This is the 
case for the so-called <r(600) scalar, often used in nuclear potentials, which has a width 
of about 600 MeV. 

2. The fact that a (980) and / (980) decay also into KK despite the highly reduced phase 
space (the decay is kinematically forbidden over most of the widths) is an argument 
against their being qq states. Indeed, four-quark states would much more readily fall 
apart to qsqs = KK than would qq scalars. In principle, the a and / could be 
"molecular", lightly bound KK threshold states, in analogy with the X(3872), which 
may be a D*D threshold state [8]. For states of similar size, the kinetic energy in 
the D*D system is four times smaller than that of the KK system. On the other 
hand, roughly the same meson-meson potentials are generated by couplings of the 
light quarks. Therefore, binding KK to form a (980) and / (980) seems unlikely. The 
features 1 and 2 above, which are particularly puzzling in a qq picture, can conceivably 
be resolved if one notes the special role of t'Hooft's anomaly-induced uuddss six-quark 
coupling [9]. 

3. Further indirect support for the four-quark picture comes from the suggestion [10] that 
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in collision or decay processes with few initial quarks, qq meson production should ex- 
ceed considerably that of more complex baryonic and exotic four-quark states. Com- 
parison of a (980) and / (980) with bonafide qq states such as p(770) mesons in e + e~ 
or pre collisions and in B decays suggests that the former are significantly suppressed, 
again supporting the four-quark hypothesis. If the initial state has many quarks and, 
in particular, many ss pairs, as is the case at the Relativistic Heavy Ion Collider, 
then the suppression of qqss production is expected to be weaker. This may be easier 
to test for / (980) than for a (980), whose identification requires good photon recon- 
struction. As further example, we note that 11% of the decay Df — > K + K~ir is due 
to f -> K+K- [5]. 

If a (980) is indeed a four-quark state, then B will be smaller than the value predicted 
utilizing ao(980) dominance and assuming it is a qq state, Eq. (1). If a search for r" — > ix~r\v T 
that is sensitive to a branching fraction of order 10~ 5 fails to detect a ~ 50 MeV-wide 
peak around 980 MeV in the ?77r _ invariant mass spectrum, this would constitute a fourth 
argument in support of the four-quark view. Conversely, observation of a clear peak would 
strongly suggest that ao(980) is in fact a regular ud state, as early arguments by Bramon 
and Masso have suggested [11]. 

Next, we present some general arguments regarding the expected scalar (L = 0) contri- 
bution Bs to the branching fraction B, assuming that it is dominated by the exchange of the 
do (980), which is taken to be a ud state. Key to its small magnitude is the operator equation 
expressing the fact that the weak vector current is conserved up to small electromagnetic 
and m d — m u mass difference effects: 

WV^x) = (m d - m u )u(x)d(x) + eA^Jx^x). (27) 

The contribution of the electromagnetic interaction term to r — > txt]v t is related to r — > 
TTfjUr'-y, but given the difficulty in observing r — > nr\v T , there is little hope that the decay 
involving an additional photon in the final state can be studied in the near future. The 
corresponding one-loop electromagnetic corrections are supressed by a/n ~ 1/500. The 
first term of Eq. (27) is ~ (m d — m u )/m h ~ 1/200 for (m d — m u ) ~ 4 MeV and a typical 
hadronic mass of m h ~ 0.8 GeV, hence we focus on this term in what follows. The matrix 
element (0|V M V^|/i) (with h = r/n or h = Oq (980), if (980) dominance holds) of the 
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operator equation (27) then yields 

Q» (OlV^h) = Q 2 f (s) = (m d - m u ) (0\S~\h) , (28) 

where S~ is the scalar current u(x)d(x) , and Q 2 = s = ml is the squared mass of the 
hadronic system. The left-hand side of Eq. (28) yields the middle expression by using Eqs. (2) 
and (3). Thus, computing Bl=o, the L = contribution to B, reduces to estimating the low- 
energy hadronic parameter (OlS 1- ^). A first-principles, unquenched lattice QCD calculation 
is lacking at present, but recent progress in dealing with light quarks/pseudoscalars may 
soon make it feasible [12]. The calculation is circumvented in the chiral perturbation theory 
approach, which uses effective Lagrangians (including isospin violation) and couplings fitted 
together to known low-energy processes and extrapolated to the r decay of interest. The fact 
that as many as three calculations of this type yielded the same result (Eq. (1)) indicates 
that this is a well-defined framework, but does not test its reliability. 

Here we present a simpler quark model-motivated estimate. Unlike the A ~ V and S ~ P 
chiral symmetry-motivated relation, we relate the axial and scalar matrix elements, since 
both pertain to P-wave (ai(1260) and a (980)) rather than S-wave (p and ix) qq states. We 
assume that (980) dominates the r~ — > ix~r\v T decay and that the decay r — > 7r~7r + 7T~i/ r 
is dominated by the a]" (1260). Defining the matrix elements 

v = (0\S-\ao (980)) , 

a = (0|A|ar(1260)i>, (29) 
where % is a helicity state index, we expect 

B L = Q l.3 V -J md - mu )\ (30) 



B(t -> (1260) i/ r ) ' a 2 V«V ( i26o) 
where the 1.3 enhancement is due to the larger phase space for the decay into the lighter 

(980). The couplings of the local scalar and axial currents to the two 3 Po an d 3 Pi ud states 

of similar mass are expected to be roughly equal, namely, a ~ v. Indeed, these couplings 

are fixed by quark-model wave functions which, apart from relatively small L ■ S effects, are 

the P-wave ground states of the same Hamiltonian. From Eq.(30) we find 

B L=0 ~ 1 x 10~ 5 , (31) 

similar to the contribution of the p~ and, within our crude approximations, consistent with 
the chiral-perturbation-theory estimates. We note that Eq. (31) may require an additional 
suppression factor of up to ~ 3, due to the three helicity states available to the a± (1260)). 
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V. TEST FOR NEW WEAK INTERACTIONS 



The general Lorentz- invariant "current x current" weak interactions could include, in 
addition to (V — A) ■ (V — A), products of scalar (S), pseudoscalar (P) and tensor (T) 
"currents". Exchanging new, heavy elementary particles cannot generate the non-minimal 
T part, hence we focus on the S and P parts. Experimentally, the amplitudes of the V ■ V, 
V ■ A, and A ■ A current products can be compared with those of S ■ S, P ■ P or S ■ P 
terms in nuclear beta decays involving both u — > d and e — > v e weak transitions [13]. It is 
convenient to parameterize the corrections to the Standard-Model currents using the same 
weak coupling g^, attributing the smallness of the S ■ S, S ■ P and P ■ P terms to heavy 
(pseudo-) scalar mesons with masses mp,m s ^> m w . A positive result implying ms,mp 
masses smaller than O(TeV) would motivate searching for such particles at the upcoming 
LHC. 

A stringent limit on the pseudoscalar mass M P comes from its contribution of g^y/Mp to 
the amplitude A(n~ — > e~v e ). The branching fraction for this decay, (1.230±0.004) x 1CT 4 , is 
in agreement with the expectation of the standard electroweak model, where its small value 
is due to the m e /m M ~ 1/200 helicity suppression of the V — A amplitude. We therefore use 
the error of this result to obtain an approximate limit on the pseudoscalar contribution, 

( Mk) 2 < 0.004 x 10' 4 — ~ 3 x 10- 6 . (32) 
\M P ) 200 v ; 

In order for measurements using unsuppressed nuclear beta decays to compete with this limit, 
a precision of about 3 x 10~ 6 is needed. Similarly, the decays K~ — > e~u e and B~ — > e~u e 
yield stringent bounds on pseudoscalar couplings involving second- and third-generation 
quarks [14]. We note that direct production of a pseudoscalar with mass Mp > 10 3 Mv^ is 
far beyond the reach of the LHC. 

The case of the scalar part is different. Current limits from high-precision nuclear beta- 
decay experiments will continue to be unchallenged by accelerator-based experiments, un- 
til an eventual 5-factory limit on or observation of the decay t~ — > ix~r\v T , whose small 
Standard-Model branching fraction makes it sensitive to new scalar interactions. In a nu- 
clear beta decay, the distribution of the angle beween the neutrino and the lepton is 

777 

W{&) = l + b-^ + af3 e cos(9), (33) 
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where m e , E e , and (3 e are, respectively, the electron mass, energy, and velocity. The beautiful 
new experiments using traps to also measure with high precision the recoil velocity of the 
daughter nucleus have observed b = -0.0027 ± 0.0029 [15], a = 0.9981 [13]. The 
deviation of a from the V — A prediction a = 1 leads to the (so far relatively weak) bound 
on the scalar mass 

~ (0.004)" 1 / 4 ~ 4. (34) 

M w 

A tighter bound of (M s /M w ) > 6 — 7 is expected from improved measurements of a. Once 
the lower part of the beta spectrum is more precisely measured, the overall normalization 
of the rate will yield a more sensitive bound of Ms > 15M^ by utilizing interference of the 
S and V — A amplitudes [16]. 

In passing, we note that standard beta decay experiments such as KATERIN [17], which 
will measure the electron-neutrino mass (or rather m Vl ) down to 0.4 eV, will have very 
high statistics of ~ 10 11 events. Still, beta spectra with or without recoiling atoms are 
also affected by radiative and hadronic effects, and precise calculations of the latter will be 
required if the experimental precision is to yield strong limits on non-standard couplings. 

A scalar ud weak current contributes to G-parity-violating second-class-current transi- 
tions, such as t~ — > it~riv T , provided that it couples to r~ and u T . As discussed above, 
the present experimental upper bound on the branching ratio for this mode is an order of 
magnitude greater than the estimated Standard-Model contribution ~ 10~ 5 , which is at the 
level that may be detected by the BABAR and Belle experiments. Since interference between 
a non-standard contribution and the small V — A amplitude will not contribute much, a 
limit of the branching fraction at the level of 3 x 10~ 5 would imply 



M s 



> (3 x 10- 5 )- 1/4 ~ 12, (35) 



M w 

comparable to the expected future bounds from beta decay experiments. 

Unlike the universal gauged weak interactions, the scalar couplings could discriminate 
between different lepton generations. Thus, the S particle could be "first-generation ori- 
ented", coupling to the u and d quarks and the e and v e leptons but not to r or u T . In 
such a case, it will affect the beta decays but not the r" — > ix~r\v T decays. Conversely, S 
particles may couple more strongly to the third-generation tu t vertex than to eu e . Thus, 
a-priori, the limit from nuclear beta decays and the one from the r~ — > ix~r\v T decay are 
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complementary and, furthermore, observation of S'-coupling effects in one mode and not the 
other would indicate non-universality. 

On the particle theory side, many lines of argument [18] suggest that new physics, par- 
ticularly novel weak couplings different from standard V — A, will most strongly manifest 
in higher generations. This would enhance S effects in the r decays relative to the first- 
generation beta decays. More generally, Ms is unlikely to be much smaller than Mp, for 
which the very strict bound above applies, unless M s is protected by SU(2) L , namely, S 
couples to the Z°. In that case, the S + , S~, and S* form an SU(2) L triplet, helping produce 
S particles at the LHC via an intermediate Z° or W ± . Otherwise, production of S + S~ pairs 
is smaller by {ocemi l( ^Weak) 2 ~ 10~ 2 . In general, if we have left-right symmetry at relatively 
low scales [19] the stringent limits on M P push M s to very high values. 

VI. CONCLUSIONS 

We have considered the S- and P-wave contributions to the branching fraction of the 
decay t~ — > ix~r\v T . We find the P-wave contribution, which is more robustly calculated, to 
be 3.6 x 1CT 6 , and the S'-wave one to be around 1 x 1CT 5 , both in agreement with previous 
calculations. Given the capability of experiments at the B factories to measure or set a 
limit on the branching fraction £>(t~ — > 7r _ ^z/ T ) at the 10~ 5 level, it is interesting to note 
the implications of the possible experimental results: 

• A "minimal" result of B ~ (0.2 — 0.4) x 10~ 5 with the 7r~r] invariant mass around the 
p~ peak, which may be hard to extract experimentally, involves no new surprises. 

• A larger value of B, in the range (1 — 1.5) x 10~ 5 , consistent with the chiral perturbation 
theory calculations and with our quark-model estimate, would strongly suggest that 
% (980) dominates the S-wave part of the decay. In this case, a narrow invariant-mass 
peak around 980 MeV should be seen. This would strongly suggest that the a (980) 
is a ud scalar meson after all. 

• A somewhat larger value, B > (2 — 3) x 10~ 5 with scalar-meson dominance, may 
indicate novel scalar components in the weak interactions. 
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